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We revisit the mean-field treatment of photoassociation and Feshbach resonances in a Bose- 
. . . Einstein condensate previously used by various authors. Generalizing the Cherny and Shanenko 

' approach ( Phys. Rev. E 62, 1646-59 (2000) ) where the finite size of the potentials is explicitly 

, introduced, we develop a two-channel model for a mixed atomic-molecular condensate. Besides the 

. individual dynamics of the condensed and non-condensed atoms, the model also takes into account 

their pair dynamics by means of pair wave functions. We show that the resulting set of coupled 
I equations can be reduced to the usual coupled Gross-Pitaevskii equations when the time scale of the 

. pair dynamics is short compared to that of the individual dynamics. Such time scales are discussed 

' in the case of typical photoassociation experiments with cw lasers. We show that the individual 

dynamics plays a minor role, demonstrating the validity of the rates predicted by the usual models 
describing photoassociation in a nondegenerate gas. 
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I. INTRODUCTION 



j3 ' Thepossibility of transforming an atomic condensate into a molecular condensate is presently a major challenge 
[E Hi IE Hi IE 1^ • Several routes are considered to couple a condensate of free atoms with a condensate of molecules 



which are all in the same vibrational state of a molecular electronic potential. A Feshbach resonance in the electronic 
ground state can be swept by a time-dependent magnetic field and recent experiments [ElEl have observed oscillations 
Q I in the number of atoms in the condensate. Alternatively, the photoassociation process, which can be considered as 
Q [ an optically induced Feshbach resonance, creates a molecular condensate in an excited electronic state. In the latter 
case, a stabilization process must be introduced to avoid destruction of this condensate by spontaneous emission as 
■ was observed recently A two-pulse STIRAP (stimulated Raman adiabatic passage) scheme has been theoretically 
^ , discussed |8„ in view of transferring the population to bound levels in the ground electronic state. 

in ; 

lO [ In a nondegenerate gas, ultracold molecules are formed through combination of a photoassociation step with a 
' stabilization step using spontaneous emission {11. 12, 13). In both cases the efficiency is controlled by the dipole 
I transition moments, which depend markedly upon details of the electronic potential curves : the search for efficient 
- mechanisms relies upon accurate spectroscopic data [T^ll5l | and it was shown that whereas photoassociation is efficient 
at large internuclear distances, the stabilization process is governed by the probability of finding the two atoms at 
"j--; ■ intermediate distances |l5l Il6| |. Making ultracold molecules thus involves an interplay between long-range and short- 
ed ■ range dynamics. Up to now, most experiments have used cw lasers. The use of chirped pulses, ie laser pulses with 
^ ' a time-dependent frequency, could increase the formation rates From the theoretical point of view, sweeping a 
I ■ time-dependent magnetic field through a Feshbach resonance is equivalent to photoassociating with a chirped pulse. 
' O ' In a condensate, most theoretical treatments at present rely upon coupled Gross-Pitaevskii equations, where the 
^ . dynamics in the atomic and the molecular condensates as well as the coupling between them, are described by 
mean- field effective potentials [E IE • Details of the potentials are omitted owing to a delta function approximation 
(contact potential). The knowledge of the molecular potential curves and dipole transition moments is necessary 
only to determine scattering lengths, binding energies and to accurately compute the coupling parameter between 
atomic and molecular phases. The success of such calculations relies on the very short time scale of the microscopic 
; ^ ' dynamics compared to the condensate dynamics. 



The validity of the (one-body) mean-field approximation has been recently questioned by several authors, particu- 
larly in the case of a time-dependent coupling term, and models using Hartree-Fock-Bogoliubov equations of motion 
have been proposed T^. However, such calculations rely on the delta function approximation for the potentials and 
coupling terms. The ultraviolet divergence caused by this approximation is solved by a renormalization procedure, 
as discussed in detail by Kokkelmans et al. ,20i|. Recently, Cherny and Shanenko 0, 12^ have shown that in the 
description of the dynamics of an atomic condensate, issues associated with the contact potentials can be avoided 
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FIG. 1: Schematic representations of two coupled channels, (a) coupling with a photoassociation laser of frequency uj/2tt: the 
asymptotic separation between the dressed potentials is given by a detuning A = Et + Ec — 2Ea — huj. (b) Feshbach resonance: 
a magnetic field is adjusted to couple the two channels by hyperfine interaction. The potentials are separated asymptotically 
by a detuning A — Et + Ec — 2Ea ■ 

by considering the exact potential and pair wave functions having the correct nodal structure at short interatomic 
distance. 

The aim of the present work is to revisit the problem of coupled atomic and molecular condensates in the framework 
of a Cherny treatment. For the sake of clarity, we shall consider only a two-channel model for Feshbach resonance 
or photoassociation. The paper is organised as follows. We first present a three-field model describing a two-channel 
coupling in a Bose system. We then derive one-body and two-body mean- field equations, and show how the pair wave 
functions can be eliminated and lead to effective one-body mean-field equations, without using contact potentials. In 
the last section, we interpret the one-colour photoassociation of a BEC in terms of one-body and two-body modes. 
Definitions of the pair wave functions in an inhomogeneous system are given in the Appendix. 



We consider three kinds of atoms in the atomic-molecular system: the ground-state atoms (a), colliding in the 
lower channel and the atoms (& and c) of the bound pairs in the upper channel. As described in Figs^ the free 
atoms a, interact through the potential Uaaif), while the molecules vibrate in the potential Ubdr). In the case of 
photoassociation, the latter potential corresponds to an excited electronic state of the molecule. For each species, we 
define a quantum field: ipa, ipb^ and ^pc■ 

This description of the system is a phenomenological starting point. It would be rigorous if the states a, 6, and c 
corresponded to well-defined atomic internal states, which is true only at large interatomic distances. Indeed, when 
two atoms come close to each other, their internal states change because their electronic clouds overlap. This means 
that, for instance, the ket -^'^i{^)ipi{y)\0) is related to the "molecular" ket -^(|x, y) -I- |y,x)) (g) where \4>aa) 

is the molecular electronic ket depending on the distance |x — y|. 

Bearing in mind this phenomenological aspect, we can write the many-body Hamiltonian of the system in terms of 
Tpa, tpb, and tpc- 



II. THREE-FIELD MODEL 



A. Many-body Hamiltonian 






(1) 
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where m is the mass of the atoms and Ea, Ei,, Ec are the internal energies of isolated atoms (see Fig^. V is 
the potential trapping the atoms, Uaa {Ubc) is the interaction potential between atoms of the lower (upper) channel, 
and Hint is the matrix element coupling the two channels (which can be time-dependent). No specific approximation 
is made regarding these potentials, so that the Hamiltonian is built on microscopic grounds. In this respect, our 
approach is very close to that of 23], the major difference being that we consider the two channels explicitly. Note 
that there are no terms involving the potentials Uab, Uac, Ubt, etc: since the atoms b and c are bound, we neglect 
their collisions with other atoms. Moreover, three-body potentials, as well as spontaneous emission (in the case of 
photoassociation), are not taken into account in this Hamiltonian. 



B. Dynamics of a two-atom system 

One can easily derive the usual dynamics of a two-atom system from the many-body Hamiltonian. For instance, 
consider the two-body state: 

\n)= Jl rf^KdV (<i>^(x,y,i)^iM|lM+ ci>A,(x,y,t)V4(x)^t(y))|o) 

where and $m are the two components of the two-body wave function for the lower and upper channels. As we 
will consider continuum states in the lower channel and bound states in the upper channel, A stands for "atomic" 
and M for "molecular". This state evolves according to the Schrodinger equation ift^^j^ = H\fl). Using the canonical 

commutation relations , {y)] = 6ijS'^{x—y) and ['0i(x), ?/ij(y)] = 0, one finds a set of equations describing a 

general two-channel coupling: 

= (^- ^ ^^2!^'^ + + ^(y) + '^-(x-y) + 2E}j <^>A + i/4t(x,y)$M (2) 

= (- +^^) + y(x) + V{y) + UU^-y) + E^ + e}\ $m 
at \ 2m / 

+ H,nt{^,y)^A (3) 

Usually, one can separate the motion of the centre of mass R = ^^^^ and only the relative coordinate r = x — y is 
considered: 

^ ('_^ + ^^^(r)+2i?„')$^(r,t) + F4,(r)$M(r, t) (4) 



^f^ d<^Mjr,t) ^ ^ Ubc{v)+Eb + E, \ <i>M(r,<) + H,nt{v)<^A{v,t) (5) 



dt 

dt \ m 

C. Dynamics of a many-atom system 

In the more general case of a many-body state, the full dynamics is given by the equations of motion for the field 
operators in the Heisenberg picture. These equations are obtained from the Heisenberg equations — 
using the canonical commutation relations. In the case of photoassociation with a continuous laser, we may actually 
rotate the field operators and use the rotating field approximation in order to remove the oscillatory time dependence 
of Hint ■ In either case, we find: 

+V2 /dV ^liy)H:nt{^,y)M^)My) (6) 



^^dM^ ^ L^^^^j^^vi.) + v(y) + uu^-y) + nA]M^)My) 

at \ 2m I 

+ -Hi«t(x,y) (7) 
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where underlining has been used as a notational convenience for symmetrizing certain quantities; for instance, 
A(x, y) actuaUy means ^(^(x, y) + A{y, x)). In Eq. iQ, we introduced the "detuning" A between the two asymptotic 
curves (see Fig^l and we neglected many-body terms corresponding to collisions or coupling with atoms external to 
the pair considered. Keeping these terms would be inconsistent with the fact that we neglected other potentials such 

as Uab, Uac- 



D. Coupling with a single bound-state 

Usually, the interaction is tuned to couple the atoms in the ground-state channel to a single stationary bound state 
in the upper channel: a rovibrational level of the potential Ubc with binding energy Etound- The relative motion of 
the bound atoms is then described by the stationary wave function ipM satisfying: 

+ Ubcir)]ipM{r) = ~EboundVM{r) (8) 

V TO / 

Here the zero of energy is set to Eb + Ec- To achieve population of this single level, one must remain in the 
perturbative regime, where the typical intensity of the coupling Hi„t remains smaller than the energy splittings 
between the molecular levels. This condition has been discussed in detail by Vatasescu et al in the case 

of photoassociation in a trap of cold alkali atoms, considering various c.w laser intensities and detunings A and 
comparing the Rabi period to the classical vibrational period of the molecular level. 

When ifiM is indeed the only bound level affected by the coupling, we can approximate the operator ^f,(x)?/;c(y) 

by its projection -0m(^^)</5a/(x— y) onto this bound level, where the "time-dependent coefficient" iPm(R-) — 
J d'^r (p*j^i{r)iipb(R- + §)V'c(R- — §) defines a molecular field. The description of the system can then be reduced to 
two fields '0a and i/'a/j satisfying: 

+ V2 Jd^ W*{^,y) 4'l{y)^M{^) (9) 

^n'J^ - (-^ + 2nx) + M)^M(x) 

A^(x+^,x-^) -^°(-+^^°(-^) (10) 

where W^(x, y) = i?mt(x, y)(/3^^(x— y) is the interaction kernel and hS = HA — Ebound is the energy detuning 
between the two levels (see Figs. considering the dressed picture for photoassociation. 



III. EFFECTIVE MEAN-FIELD THEORY 



A. Purely atomic system 

Let us first consider a purely atomic system (thus ignoring the terms involving any field V'f) or ipc). The equation 
of motion © now simply reads: 

^^dj^ - (-^ + + /^V cUx-y)^l(y)V3.(y)) U^) (n) 

In the condensed phase, one usually assumes a non-zero expectation value 5'o(x) = {ipa{'x.)) for the field operator 
-00 (x), corresponding to a macroscopically occupied state (see the Appendix). What is sometimes referred to as 
"naive mean-field" consists in replacing the field operators by their averages directly in the equation of motion l|ll|l . 
thus neglecting the quantum fluctuations 6a = i/'a ^ ^o- This leads to a non- linear Schrodinger equation with a 
coupling constant of the form Jd^ Uaa{'^)- This coupling constant is well defined for weak-coupling interactions 
(for instance, interatomic potentials satisfying the Born approximation Jd'hc Uaa{'^) ~ ^^^^r~ ■'^h^re a is the s-wave 
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scattering length). It is not the case however for the interactions considered here. The interatomic potential exhibits 
a strong repulsive hard-core which leads to a divergent coupling constant. 

The usual remedy is to replace the real interaction by an effective one |23|, generally a contact potential ?7Qa(r) — 
4-Kah gs^j.^ (which does satisfy the Born approximation) having the same scattering length a. The resulting equation, 
known as the Gross-Pitaevskii equation (26> . i27i | , forms an "effective mean- field" theory, in which only the large-scale 
effects of the interaction are retained: 

.^d*o(x) / fi^Vl . Airah^ . , ..^ ^ , , , ^, 



This method can be rigorously justified by considering the effect of two-body correlations |2lLl28l |. An intuitive way 
of taking these correlations into account is to use pair wave functions. Just as most atoms are individually described 
by the same macroscopic wave function ^o(x), most pairs of atoms are described by the same macroscopic pair wave 
function <i>oo(x, y), which turns out to be simply the anomalous average -^(i/'a(x)?/;a(y)) (See Eq. (|A18|) in the Ap- 
pendix). At large distances, the atoms are decorrelated, and the pair wave function is just a product of two macroscopic 
one-body wave functions (^/ia(x)Vio(y)) « {tpaix)) {ipaiy)} ■ So we may write $00 (x, y) = ^^'o(y)*o(y)</'oo(x, y) with 
lim|x_y|^oo 'Poo(x, y) = 1. The function (poo may be seen as a reduced pair wave function describing the correlations 
at short distances (see FigE)). To some approximations (see assumptions HI and H2 in the Appendix) justified by 
the low density and large extent of the condensate, ipoo simply satisfies the scattering equation for two atoms in free 
space. Written with the centre-of-mass and relative coordinates R and r, this equation reads: 

^^cWR^ ^ (-£v|,,-|^V? + f/..(r)) ^oo(R,r,t) (13) 



dt 

At equilibrium, it becomes a stationary state tp^-^'> satisfying : 

E^(^\R, r) = ( - ^ + [/,,(r)) ^(^)(R,r) (14) 

When the condensate pairs are in their ground state (which is presumably the case when the condensate is not 
excited), (poo is therefore the lowest energy state satisfying lim|r|^oo '^^^■'(R, r) ~ 1, ie the stationary scattering state 
at zero energy ip^^^ |42]. This state can be expressed formally by: 



\r) = l- J rf^R'dV Uaa{r') G{-R, r, R', r') - 1 - Jd\' Uaa{v') g{v, r') (15) 

where G is the Green's function of the operator — ^V|j — Tn'^'r + C^aa(r) and g the Green's function of the operator 
— — -|- C/aa(r). It is known from collision theory that this state is related to the s-wave scattering length a: 



d\Uaair)^^'Hr)^^^g (16) 



When taking the expectation value of the equation of motion (|lll) . one finds the term {'ipl{y)ipaiy)i^a{^)} , which can 
be approximated by (^I(y))('?/'n, (y)'0a(x)) = |^'o(y)pvE'o(x)iy9oo(x, y) when neglecting the collisions of non-condensed 
atoms with condensed atoms {pa^ ^ 1) |2ll |. If we assume again that the averaged field ^'o is nearly uniform on the 
scale of pqo, the interaction term becomes at equilibrium: 

(^J dV L/aa(x-y) ^(°)(x-y)^ |*o(x)|2vl/o(x) = l!^|vl/o(x)|2M.o(x) (17) 

where we used H16|) . Thus we retrieve the mean- field term in the Gross-Pitavskii equation l|12() : it is in fact 
regularised by the stationary two-body correlations described by p'-^K 

B. Atomic and molecular system 

We now apply these ideas to the atomic-molecular system, with the aim of deriving an effective mean-field theory. 
A first approach ITsL l30l| would be to start from the two-field description l|^. (|10|) and replace the potentials by 
effective (or renormalised) interactions: C/aa(r) — 3'5'^(r) with g = '^'"^ ° and W{v) — 'w5'^{v), with 

w = ^/2{p>M\H,ntW^''^) = V2 jd\ ip*M{r)H,„tir)p'-°\r) (18) 
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Replacing the quantum fields by their averages {tpa) = and (M ^m) = ^'m, we obtain a set of coupled Gross- 
Pitaevskii equations which has been extensively studied 0, IsL llsL IsoL Isil l3^ : 



,.d*o(x) 




h ; = 

dt 


\ 2m 




( f^'^l 


dt ^ 


\ 4m 



F(x)+5|^'o(x)|^j^'o(x) + «;*^'*(x)*m(x) (19) 
2F(x) + *Af(x) + i«;*g(x) (20) 

Let us now follow the approach of ,21j. We now consider the following one-body and two-body fields: 

• the atomic condensate mode ^E'o 

• the atomic non-condensed modes "ili 

• the molecular condensed mode $a/ = {fpb'4'c) (we will neglect non-condensed molecular modes) 

• the pair wave function $00 for two condensed atoms 

• the pair wave function $oi for & condensed atom and a non-condensed atom in the mode ^'^ (we will neglect 
the pairs of non-condensed atoms) 

We first take the expectation value of Eqs. ©,0, expressing the atomic correlations by means of pair wave 
functions: 

^f^^^M. ^ ( - ^ + ^(x))^o(x) +V2jd^ ( Uaa{x-y) vI':(y)<i>o.(x,y) 

i 

+ 7l4,(x,y)vE'5(y)$M(x,y)) (21) 



^n'-^^ = (-^^!^^±^ + nx) + ny) + t/.(x-y) + ^A)<I>M(x,y) 
at V 2m / 

-f iI™t(x,y)$oo(x,y) (22) 

The dynamics of 5*^, $00 and $oi is given in the Appendix. Note that Eqs. (|A34p and l|22() giving the evolution of 
$00 and are analogous to Eqs. Q,©. They describe indeed the coupled dynamics for the pairs of condensed 
atoms, giving rise to Rabi oscillations on a time scale 

Tpairs — • (^^) 

n-int 

To simplify these general equations, we make use of the reduced pair wave functions (^00 and (^oi defined by Eqs. 
(|A11() . HA12() . with the assumptions Hi and H2. Furthermore, we resort to a perturbative approach and assume that 
the atoms are coupled to a single stationary bound state (^sm, satisfying Eq. (|HJ|. In this case, $m(x, y) reduces to its 
projection ^'m(^^^)va/(x— y) onto ipM- The "time-dependent coefficient" m is the centre-of-mass wave function 
for the molecules in this bound state. 

These assumptions lead to a closed set of equations: 

ifi-^ = { - + V + gool^of + J22go,\^^\^)^o + 3m*m1'S (24) 



dt V 4m m 



+ v^H.-..(R.r)^„|r)( ^^^^;,';g|^_,^ ) (27) 
aif^ . (-^!Zl-^!p + c..„),„.(K,,) (28) 
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where the one-body modes are coupled to the two-body modes through the couphng factors: 

goo(x) = J(fy t/aa(x-y) y)oo(x,y) (29) 

5oj(x) = Jd^ Uaa{x-y) Lpoi{x,y) (30) 

5m(R) = V2y"A (^;,(r)i/„t(R,r)(^oo(R,r) (31) 

5;^(x) = \/2 ^;^(x-y)g,„t(x,y) (32) 

Note that the non-condensed modes as well as the pairs tpoi involving a non-condensed atom arc not directly 
affected by the molecular condensate. Thus, our equations cannot present any "rogue photodissociation" |30l l33| . 
enabling molecules to dissociate towards non-condensed modes. The rogue photodissociation may be found either in 
the fluctuations of the molecular field or in the pairs ipij involving two non-condensed atoms, which we have both 
neglected. However, rogue photodissociation has been treated thus far without the molecular fluctuations, and using 
the anomalous average (aka_k) to describe non-condensed pairs. In our own approach, this average is rather related 
to the condensed pairs at short distances (see Eq. (lAlSp ). We plan to clarify this point in a future work. In the rest of 
this paper, we will assume that we are in situations where the rogue photodissociation does not play an important role. 

In order to decouple the one-body dynamics from the two-body dynamics, the pair wave functions must be quasi- 
stationary and follow adiabatically the one-body wave functions. In other words, the characteristic time of evolution 
for the pairs Tpairs (see Eq. must be very short compared with the characteristic time of evolution Tcond for the 

condensates: 

« 1 (33) 

^ cond 

In that case, we can solve Eqs. (|27|l . H28f) formally using the Green's functions introduced in Eq. (|15|1 : 
-J d^R'dSr' G(R',r',R,r)V2i/,;,i(R',r')¥'A/(r') 



^m(RO 
*o(R' + t)*o(R- 



■J 



„^Af(R) 

(R) 



¥>o.(R,r) = ^("Hr) 

We can then eliminate </3oo and iy9oi in Eqs. |(53J , ((231 , (EEl , and using the properties lfT ^ . (fT?))l of (/5^°\ we find: 

= I -^+y + 5(|*o|' + ^2|*,|2) I ^ y;*^*^^, (34) 

..^^ (-^ + V + 2,|^oP)v.. (35) 
= (^-^ + 2V + h6 + E,,if'^^SM + Iw^l (36) 

where Eseif = / K-^^I^^^'I-p )I p(E}dE is the self-energy of the molecules: it is an energy shift caused by the 
interaction with the open channel 0, Hi El El . We may include this energy shift in the detuning 6. Note that 
Eq. (|35|l is conservative, so that the non-condensed modes may be omitted if they are initially negligible. Finally, we 
retrieve the coupled Gross-Pitaevskii equations (|19I2U|) . 

One can see in Eqs. that the typical time Tco7id for the Rabi oscillations between and \E'm is given 

by h/{wy/p), where p is the typical density of the system while the interaction w has been defined in Eq. p8(l . The 
condition (|33|l is thus equivalent to: 

MvMlip'^"^) « 1 (37) 
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IntGiisity^ 


Eh ^ 

Abound 


T 

^ pairs 


cond 


^ spont 




'^'■"Na \7] 


100 W/cm^ 


49 cm,-^ 


8.4 X 10-'"s 


1.1 X lO-'-s 


8.4 X lO^^'s 


7.4 X 10"* 


1 kW/cm^ 


49 cm~^ 


2.6 X 10"'"s 


3.5 X 10"'s 


1 kW/cm-' 


4 cm^^ 


2.6 X 10"'"s 


3.5 X 10^'s 


4.5 X lO"'' 


"'"Rb \8] 








4.8 X lO^'^s 


1.4 X IQ-'^s 


~ 3.10"* 



TABLE I: Some typical time scales for the photoassociation of a condensate. 



Note that this condition does not depend on the intensity of the coupling Hintj but keep in mind that it has been 
derived in a perturbative way. More generally, mean-field equations such as H34|) . H3t)|l should be valid as long as the 
time scale Tpairs, associated with the two-body dynamics, is short compared with the time scale Tcond of the one-body 
dynamics. In other cases, one might need to use the more general Eqs. and treat the two-body dynamics 

explicitly. 

C. Comparison with other models 

These general equations can be compared with the models developed by Holland et al. 19] and Kohler et al. [23l |. 
In the Hartree-Fock-Bogoliubov model of Holland et al., anomalous and normal correlation functions are used to go 
beyond the coupled Gross-Pitaevskii equations: when these correlation functions are negligible, their model also leads 
to the coupled Gross-Pitaevskii equations. We therefore believe that the anomalous correlation is in fact related at 
large distances to the fluctuations around the stationary perturbed pair wave function $00. It is not true, however, 
at short distances, since their model is built with contact interactions. This is why a renormalization procedure is 
needed to recover the physics from short distances. 

In the model of Kohler et al, the system is described by a single non-linear non-Markovian Schrodinger equation 
which is not obviously related to the coupled Gross-Pitaevskii equations. However, the cumulant approach at the 
basis of their model should be mathematically very close to our pair wave function approach. The main differences 
are that we factorize the correlations and consider the two channels explicitly: what they note \E'(x)'I'(y) -|- $(x, y) is 
a condensed pair wave function with two implicit components corresponding to the pair wave functions $00 and $ m ■ 
As our model has an explicit connection to the coupled Gross-Pitaevskii equations, we believe it is also the case for 
the model of Kohler. Our model therefore clarifies the relation between the existing models and gives a simple image 
in terms of one-body and two-body wave functions. 

IV. APPLICATION TO PHOTOASSOCIATION IN A BEC 

The condition H37|) shows that the one-body mean-field approach holds for sufficiently low densities, contrary to 
what has been suggested in ref . 8j . Actually, in the case of photoassociation, the condition (|S7|l is usually achieved 
in the current experiments where the condensate is always dilute enough {p ~ lO^^m"^), and the Franck-Condon 
factor ((/PAf jiyS^"^) is sufficiently low ('^ lQ~^^m^/^). We give in Table I typical values. We first have considered the 
photoassociation experiment in a Na condensate reported in 0, where the level v = 135 of Na2 0+(35' -|- 8^1/2) 
(corresponding to Ebound ^ 49 cm^^) is populated for various laser intensities. The computed radial wavefunctions 
are represented in Fig.© : although the Franck-Condon factor is favourable for photoassociation, we find for a 
density of 4 x 10^° a ratio Tpairs/Tcond < 10~^. It is not easy to increase the transition moment by three 

orders of magnitude. Photoassociating towards molecular levels closer to the dissociation limit, for instance v = 163 
(corresponding to Ebound ^ 4 cni~^ and an outer turning point located at r ^ 85 au), would increase the Franck- 
Condon factor by less than one order of magnitude. Similar conclusions on the relative time scales can be drawn for 
a Rb condensate, in the conditions described in Ref.j^. We then expect that in most realistic cases the condition 
Tpairs ^ Tcond wiU indeed be verified. 

There is another time scale, however, since the molecular state is not stable and decays by spontaneous emission. We 
can treat this spontaneous emission phenomenologically within the framework of our equations by adding a loss term 
—ih^'i']\j in the equation of the molecular field H26|) . To simplify, we will consider a homogeneous system {V = 0). 

For sufficiently high intensities. Hint 3> hj, ie Tpairs ^ Tspont 1 and the pairs have the time to oscillate between 
the two channels before spontaneous emission takes place. Eliminating the pair dynamics, we are left with a set of 
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FIG. 2: Representation of the radial wave functions in the photoassociation experiment of Ref.Q- The function in the open 
channel and the bound wave function were computed using respectively the Numerov-Cooley method fsSj and the Fourier Grid 
Method Note that in the text the total wave functions (poo{r) = Foo{r)/^/4ivr and (pM{r) = FM{r)/V47rr are considered. 



coupled Gross-Pitaevskii equations: 



= 5l*o|'*o + w*%-^M (38) 
^^d^ = hid-tj/2)^M + Iw^l (39) 



dt ' '/ / 2 



Spontaneous emission is then usually faster than the dynamics of the condensates {hj ^ w^/p), so that we can 
eliminate the molecular field adiabatically: 

We obtain a simple rate equation for the condensate density l^'op: 



dt 



K\^o\' (41) 



with a rate K — , , and an on-resonance rate: 

l+(2i5/7) 



K^ = -\U^M\mntW^'''>)\^ (42) 
7 Ti 

This is exactly the rate one can derive from a perturbative treatment of the two-atom equations Q , ((SJ [34] , assuming 
the presence of the condensate has only two effects: very low collision energies in the open channel and no need to 
symmetrize since condensed atoms are in the same state. This last effect reduces the rate by a factor of 2 relative 
to a non-condensed gas |30l l37j , and was accounted for in our model through the proper symmetrization of the pair 
wave functions (compare Eqs. (|A11|) and (jA12|) '). Thus, a cautious many-body treatment of the photoassociation in 
a BEG eventually yields the same prediction as two-atom theories. This indicates that the physics of the problem lies 
essentially in the pairs of atoms. This result has been confirmed by the photoassociation experiment in a condensate 
of sodium atoms reported by McKenzie et al 0- The laser intensities used in this experiment (from ~ 0.1 to 
1 kW.cm~^) are such that we are in the case considered here, where Tpairs ^ Tspont (see Table QJ. They found 
that the condensate is indeed locally depleted according to Eq. (|41|l . with a rate Kq proportional to the intensity, 
such that dKo/dl = 3.5 x 1Q~^^ s~^W~^ crri^ . This is in agreement with the theoretical rate given by Eq. 142|l . 
Our own calculation, using numerical computations of stationary wave functions ipM and tp'^^^ (see Figl^jl), leads to 
dKo/dl — 3.3 X 10^^^ s^^W~^cm,^, in agreement with the calculation given in 0- 
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For smaller intensities, Hint ^ ^7, and the spontaneous emission is the fastest phenomenon: the formed molecules 
decay almost instantly, so that we can eliminate them adiabatically. Setting d^M/dt = in equation (|26|l yields: 



xp,, ^ - "^^Z ; Z '^i (43) 



^/2h{5 - 17/2) 



So that we are left with: 

~dt 



(</?A/|i?int|<^00),T,2 


= ( /d^r J[(poo](r))|*oP*o (44) 



ih^lv) = _^vVoo(r) + Jiy'ooKr) (45) 
dt m 



where J[(^oo](r) = f/aa (r)¥'oo(r) — -ff*„j(r)i^ji/(r)-^^^|^|j^^|^^. The situation is the following: the formation of 
excited molecules creates a hole at short distances in the pair distribution (/Sqo i through the imaginary part of J. This 
hole is refilled with waves coming from larger distances, while the whole condensate is being depleted. Although this 
picture is quite different from the higher-intensity regime, we believe it leads essentially to the same rate of depletion, 
as we do not expect to see new behaviour emerging from lower intensities. 



V. CONCLUSION 



We have generalized the pair wave function approach, introduced by Cherny for an atomic condensate, to a 
non-homogeneous Bose system with two coupled channels, such as a mixed atomic-molecular condensate. The 
treatment of pair correlations is performed using the physical potentials both in the atomic and the molecular channel. 
We have shown that the pair wave functions regularise the mean-field terms in the equations without requirement of 
any renormalization procedure. This has enabled us to derive the coupled Gross-Pitaevskii equations on firm grounds 
and to assess their range of validity by comparing characteristic times associated to the pair dynamics and to the 
condensate dynamics. 

We have shown that in the case of photoassociation with a cw laser, at intensities allowing a perturbative 
model, the Gross-Pitacvskii description is usually verified and even leads to the rates predicted simply from the 
photoassociation probability of two colliding atoms. Nevertheless, we believe that the two-body mean-field equations 
might be necessary in other situations with a non-trivial time dependence, such as photoassociation with chirped 
laser pulses. Previous works J.9, 23, 40] have already shown the importance of long-range correlations for Feshbach 
resonance induced by a time-dependent magnetic field. Future work will clarify the possible influence of short range 
correlations in these time-dependent situations. 
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APPENDIX A: MODES AND PAIRWAVE FUNCTIONS 

In this appendix we generalize the ideas of Cherny to a non-homogeneous Bose gas. To simplify the notations, 
the field operator ipa will be noted ip. 

1. Definition of the pair wave functions 

In second quantization, the one-body and two-body density matrices of the system are defined by: 

pi(x,y) EE (^^(x)7^(y)) (Al) 
P2(x,y,y',x') ^ l(^t(x)^t(y)^(y)^(x')) (A2) 
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As hermitian matrices, they can be diagonalised in an orthonormal basis: 

Pi(x,y) = Y.n,r^{^)i^^{y) (A3) 

i 

/92(x,y,y',x') = ^m^C(x,y)0^(x',y') (A4) 

V 

This means that there are atoms in the one-body mode (or wave function) 7/;^, and there are pairs of atoms 
in the two-body mode (or pair wave function) All these wave functions are normalised to unity, but we can 
alternatively define the functions — ^/n^ij^i and — ^/rn^(f>^, which arc normalised to the number of atoms or 

pairs of atoms. Note that in a state where the total number of atoms is n, we have = n and = "^"^^"'"^ . 

This may be seen as a consequence of the relations: 



d\ (V/(x)V;(x)) = n (A5) 

(^t(x)^t(y)^(y)^(x')) ^ („-l)(^t(x)^(^')) (A6) 

We can use the one-body modes ■04 to expand the field operator: ^/'(x) — aiipi{x.). Using (jAip and (jA3|) . we find 



When there is no interaction between the atoms, the atoms are decorrelated and the pair wave functions are 



the relation (a|aj) = riiSij 



simply symmetrized products of one-body wave functions: 



0,^(x,y) = ^/>^(x)V',(y) (A7) 
ipi (xjipj (y) -I- ipj (x)'0j (y 
V2 



,,(x,y) = ^..v^v^.w. 



is the wave function for two atoms in the same mode tpi, and is the wave function for two atoms in the 

modes t/^i and t/jj. From counting arguments, we have ma = pairs in the mode (fya and = iiiUj pairs in 

the mode (pij . One can check that this pair distribution does indeed satisfy the relation (|A6(I . 

In the presence of interactions, one might expect extra terms needed to describe the correlations at short interatomic 
distances |x — y| (see Fig|31l: 

(/)„;(x,y) cx V4(x)V'»(y) + 4(x,y) (A9) 
, . ^ V'»(x)V'j(y) + V'j(x)V'4(y) , 

(?!)y(x,y) cx h0jj(x,y) (AlO) 

These correlations are expected to vanish for interatomic distances larger than, say, the range tq of the inter- 
action. T he p air wave functions are thus asymptotically decorrelated; this follows from the principle of correlation 
weakening j21j. If the interaction supports bound states, one might also expect bound-pair wave functions 0^, vanishing 
entirely for interatomic distances larger than tq. 

As the functions 4>'j^j will scale like the functions ipi and ipj, it can be useful to express the correlations with 
dimensionless functions, which we refer to as "reduced pair wave functions" : 

^■«(x,y) cx ■;/'j(x)V'j(y)'^ij(x,y) (All) 
J, I ^ '0j(x)-0j(y) +'(/'i(x)'(/'j:(y) / N 

i>^M.y) a — '/'y(x,y) (A12) 

with <y9y(x, y) — > 1 when |x— y| 3> ?'o. 



2. Pair wave functions in a condensate with U(l) symmetry breaking 

When condensation occurs, most of the atoms are in the same quantum state, say (^q. This means tiq 3> X^j^i- 
Accordingly, most of the pairs are in the same quantum state V'ooi corresponding to pairs of condensed atoms. A way 
of treating this situation, related to the Bogoliubov ansatz, is to assume a breaking of the U{V) svmmetrv|4lj: the 
system is then in a coherent state jfi) satisfying ao|ri) — ao|^^)- From the relation {d\aj) — fiiSij, we obtain ao = y/no 
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FIG. 3; Schematic representation of the macroscopic wave function for condensate pairs in a one dimensional system: 
^oo{x,y) = -^"^1 o{x)'i o{y)'fioo{x , y) ■ The reduced pair wave function (foo gives the correlation at short distance \x — y\ created 
by the scattering of the two atoms. 



and (di) — {i ^ 0), so that the field operator has a non-zero expectation value (-0) = vI/q. The deviations from 
the mean = 'I'o — V' £^re the quantum fluctuations. This implies that the number of atoms now fluctuates and the 
relation (|A6|I derived for a fixed-number state is no more valid. However, this prescription is very useful to derive 
expressions for the one-body and two-body modes in terms of the quantum fluctuations. Indeed, if we expand the 
two-body density matrix in terms of d: 



P2(x,y,y',x') = i{vI/S(x)*S(y)*o(y')*o(x') + *S(x)*S(y)(^(y')^(x')) (A13) 

+ (0~(x)^(y))*vl/o(y')vl/„(x') 

+ vE'*(x)*o(x')(%)^(y')> + *SW*o(y')(^^(y)^(x')> (Ai4) 

+ vE'*(y)vI/o(x')(^^(x)^(y')) +vI/S(y)vI'o(y')(^^(x)^(x')> 

+ vE.*(x)(0t(y)^(y')^(x')) + Vl/*(y)(^t(x)^(y')^(x')) (A15) 

+ *o(y')(^^(x)^ny)^(x')> +*0(x')(^nx)^t(y)^(y')) 

(0~t(^)^t(y)^(y')^(x'))| (A16) 



We see that it can be factorized as follows: 



where: 



P2(x, y, y', x') = $^(x, y)$,, (x', y') + ^ ^(x, y)<i>.(x', y') (A17) 



(V.(x)^(y)) ^o(x)^o(y) ^ {aCa,) ,,,,,, 

$oo(x,y) = -= = -= + 2^ — -^V'»(x)V'i(y) (A18) 

^2 V2 V2 

^ ^o(x)^.(y)y.(x)^o(y) ^ ^ i|Mi0,(x)V.(y) (A19) 



are respectively the pair wave function for two condensed atoms, and the pair wave function for a condensed atom 
and a non-condensed atom in the mode "^i. There is no simple expression for the pair wave function <i>ij of two non- 
condensed atoms, nor the bound-pair wave functions since these are contained in the term (0^(x)0'''(y)0(y')^(x')). 
If we assume that there is no interaction between the non-condensed atoms (which is the case in the usual Hartree- 
Fock-Bogoliubov approach), we may use Wick's theorem to expand this term: 

(^t(x)^t(y)^(y')^(x')> = (^t(x)^t(y))(0(y')^(x')) 

+ {0^{^)e{^)){e\Y)e{Y')) 

+ (0t(x)^(y'))(0t(y)^(x')) (A20) 
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This yields: 

*o(x)*o(y) 



$oo(x,y) = I3±ipyi+ J2 '^U^)Uy) (A21) 
v2 ,-/7rf/n v2 



1 

V2" 



cl>.(x,y) = (A23) 

$.(x,y) = (A25) 

As expected in this case, only the condensed pairs are correlated, through the anomalous correlation (6'(x)6'(y)) = 
^^_^g ^■_^g(aiaj)V'i(x)V'j(y)- This means the condensed atoms interact only with each other, while the non-condensed 
atoms are treated as an ideal gas. However, this is an approximation, since all the atoms feel the interaction. As 
we shall see, we actually need the correlations in order to regularise mean-field terms appearing in the equations of 
motion for the one and two-body modes. 



3. Evolution of one-body modes 

The evolution of the condensate mode is simply obtained by taking the average of the equation of motion for the 
field operator Doing this, we encounter the average (V'^(x)^(y)V;(z)), which can be expanded and then factorized 
in terms of pair wave functions: 

(^t(x)^(y)v;(z)) = %{^){e{y)e{z)) + ^^{y){e\^)e{7.)) 

+ *o(z)(^t(x)0(y)) + (^t(x)^(y)^(z)) (A26) 
= y2 5]vi;*(x)$o.(y,z) (A27) 

i 

We also have to express the average (V'J(y)iJ*„f (x, y)V'b(x)?/'c(y)) which can be simphfied to 
^'Q(y)i?,*„j(x, y)$M(x, y) if we assume that the atoms of the upper channel are paired in a purely coherent 
field <i>M- The equation of motion for the condensed mode is then: 

+ V2jd^ (yX^iy) Uaa{^-y) ^o^i^,y) + %{y) H■Lt{^,y) ^M{^,y?j (A28) 

We can then derive the evolution of the one-body density matrix for the non-condensed atoms pi (x, y) = 
{ijj'^{:x.)ijj{y)) - ^';5(x)^'o(y) = Z^i^^o (^)*i(y)- Using Eqs. ©, l|A17|l . (|^ . and neglecting the pair wave func- 
tions $.y for two non-condensed atoms, we get: 

^n'-^ = i^H^)^n'-M)-%i^)^n'^-{.^yr 

( fi^V^ \ 

= -^ + ny) P'i(x,y) 

+ 2^d\ <i>^;o(x,z) g*,,(y,z) j>M(z,y) - {x ^ y}* (A29) 
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In the expression above, the correlation terms $'(x, z) vanish for |y — x| vq. So when x and y are separated by 
macroscopic distances, we have: 



ih 



dt 



which can be written as: 



-^ + ny)U'r(x,y) 



+ 2 dhUaaij^Y. 



V2 



$o*(z,y) - {x^ y}* 



(A30) 



^ — -ji {x-yj 



17^0 



dt 



2rn 



V{y)] ^^{y) + V2 jdh t/,,(y-z) ^g(z)$o,(z,y) 



It is clear from this expression that the evolution of the non-condensed mode is given by: 



{x^y}* (A31) 



ih- 



dt 



2m 



V{x)^ *,(x) + V2 jd\ t/.a(x-y) ^g(y)$o,(y,x) (A32) 



as long as its wavelength is greater than tq. The interpretation of Eqs. (|A28|I and (|A32|I is straightforward. Every 
condensed atom evolving in the trap can collide with another condensed atom or with any non-condensed atoms, and 
it can also associate with another condensed atom to form a molecule. Similarly, every non-condensed atom evolving 
in the trap can collide with any condensed atom (we neglected the collisions with other non-condensed atoms). Note 
that, within the approximation of a purely coherent molecular field, the non-condensed modes are not directly affected 
by the formation of molecules. 



4. Evolution of the two-body modes 



The equation of motion for the pair wave function of condensed atoms follows from Eqs. (jA18|) and JB)): 

z;.^(x,y) = -L(V,(x)zfi|(y)) +{x-y} (ASS) 

= { - ^^^^l^^^"^ + V{^) + V{y) + Uaa{^-y) \ $oo(x,y) + i7,L(x,y)*A/(x,y) 
\ 2m / 

+ jd\ ( Ugg (z-y) a2 (z, z, x, y) + pi{z, x)g*,, (y, z)^m (y, z) + {x ^ y}) (AS4) 

with (T2(z,z,x, y) — ^(■0^'(z)V'(z)'0(x)V'(y)). Using Eq. ljA28p . one can then deduce the equation of motion for 
the reduced pair wave function (poo ■ 



^ dy)oo(x,y) 
* dt 



-(Vln^'o(x) • + Vln*o(y)-Vy) 



— (Vx+Vg) + t/aa(x-y) 



+ d^ ( Uaaiz-y) F{x,y,z) + G(x, y, z)\/2i7,„t(x, z)$a/(x, z) + {x ^ y}) (/3oo(x,y) 



/^TT* I ^ ^M(x,y) 
V 2iij„j (x, y) 

'•''l'o(x)*o(y) 



(ASS) 



with -F(x, y, z) 



and G(x, y, z) 



(i/»t(z)V-(y)) (■0t(z)) 



These many-body terms 



(V(x)V'(y)) (V'(y)) ^ ' (y'(x)^(y)) (V(x)) 

are caused by the interaction with other atoms external to the pair considered. At sufficiently low densities, we can 
make two basic assumptions about the reduced pair wave function: 
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(a) \, __,,/ (b) 



FIG. 4: Schematic representation of the effective potentials felt by (a) the condensed atoms, (6) the non-condensed atoms (for 
the first non-condensed modes), in the case of a positive scattering length. 



HI: its extent ^ is mueh smaller than the typical length scale of the condensate 

H2: it is not much influenced by remaining many-body terms like F and G. The many-body effects on the pair 
dynamics are thus only caused by the mean fields ^i, through the relations I|A11|I and l|A12(l . 

Within these assumptions, Eq. I|A35|I reduces to: 

ift =y-^^^-+^v^^ U--i^-y) ) (x, y) + V2i/,„, (x, y) ^^^^^^^^^^ (A36) 

The pair wave function $oi for a condensed atom and a non-condensed atom in the mode i is given by 
Jd^z 4'i(z)(V;'''(z)-0(x)-0(y)), following from Eq. ljA27|l . Differentiating this expression and using the assump- 
tions Hi and H2, we find for the reduced pair wave function ipoi: 

^^d^po^ ^ ^ _ |l(V2^V2)^ ^^^(^„y) j^^^^(^^y) (A37) 

Within our approximations, a pair involving a non-condensed atom is not directly affected by the coupling between 
the channels. This may indicate that one actually needs the quantum fluctuations of the molecular field for these 
modes. 

Under the assumption Hi, we can simplify the equations l|A28|l and (|A32|) for the one-body modes: 

= (-^ + ^ + <?ool*ol' + E 250.1*.!') *o + Ga/*S (A38) 



with: 



= (-^ + ^ + 2ffo.|*on*i (A39) 



500 (x) = jd'^y ;7aa(x-y) y)oo(x,y) (A40) 
go,(x) = d^y Uaaix-y)'Pot{^,y) (A41) 



Gm(x) = V2 yc;3y g4^(x,y) $M(x,y) (A42) 

These factors are proportional to scattering amplitudes. Note that the amplitudes are twice larger for pairs of 
condensed and non-condensed atoms; this comes from the fact that the atoms are in different states, so that the 
symmetrization (jAlOfl induces a factor of two relative to a pair of condensed atoms. As a consequence, the effective 
potentials felt by the condensed atoms and the non-condensed atoms are different. Fig01 gives a schematic represen- 
tation of these potentials at zero temperature equilibrium, where goo — 9oi = .9 = ° ; with a positive scattering 
length a. 
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